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We present a measurement of the total width rjjj t and mass M Vc of the r\ c meson with the 
BABAR detector at the PEP-II B-factory at SLAC. The results based on a data sample of 
88 /b" 1 are T^ t = (33.3±2.5(stat.)±0.8(sys.)) MeV/c 2 and M Vc = (2983.3±1.2±1.8) MeV/c 2 . 
An enhancement observed around 3.63 GeV/c 2 can be interpreted as an evidence for the 
t) c (2S) state. In addition a preliminary measurement of the r lepton lifetime using 30.2 /b -1 
is described. The measured r lifetime is: r T = (290.8 ± 1.5(stat.) ± 1.6(syst.)) fs. 



1 Measurement of the rj c meson mass and width 

1.1 Introduction 

The mass and the width of the charmonium ground state rj c (lS) are not well established. The 
world average d of the total width and mass is T^ t = IG.OtH MeV/c 2 and M Vc = (2979 ± 
1.5) MeV/c 2 , respectively. The rj c is expected to decay dominantly via two-gluon annihilation, 
so that r^ t ~ T^ c g . Hence the ratio of the total width of the rj c to its two-photon partial width 
is predicted by NLO perturbative 



T Tot _ 9 «s 1 + 4.8a s /?r 
~ 8o7 ' 1 - 3.4a s /vr 



(1) 



Using^l a s = 0.28 ± 0.02 and the world average^ of T 1 ^ = 7.5 ± 0.8 keV gives F^ t as 25.4 ± 
6.0 MeV/c 2 , a value higher than but still consistent with the world average. 



1.2 Fit of the mass spectrum 

The data sample consists of an integrated e + e~ luminosity of 88 fb , including runs on and 
below the T(4S) resonance. The r\ c mesons are produced by two-photon interactions. They are 
selected via the decay channel rj c — » K^K^ir^ with K® — > tt + tt~ , where in addition the e _ and 
e + are required to escape through the beam pipe. 

The mass spectrum is shown in figure ^ A large peak can be seen at the r\ c mass, a smaller 
one at the J ftp mass. J/ip are produced by e + e~ annihilation where one electron emitted a hard 
photon by initial state radiation. 

The width of the r\ c is of the same order of magnitude as the mass resolution of the detector. 
Therefore the detector resolution has to be precisely known in order to extract the r\ c width from 
the observed peak in the invariant mass spectrum. The width^of the J/tp of 87 ±5 keV is much 



smaller than the detector resolution. Hence the observed width of the J/ip peak is completely 
determined by the resolution. 

The mass spectrum is fitted in the mass range from 2.4 to 3.6 GeV/c 2 . The r\ c is represented 
by a convolution of a non-relativistic Breit-Wigner function and a Gaussian. The J/ip peak 
is fitted with a single Gaussian. The standard deviation a Vc of the Gaussian, which describes 
the mass resolution at the i] c peak, is constrained to a value 0.8 MeV/c 2 lower than that for 
the J/ip. This difference of 0.8 MeV/c 2 of the mass resolution was obtained from Monte Carlo 
events. The background is represented by the function: A ■ e ~ B - mass . The free fit parameters 
are the J/ip mass Mjm, the mass difference M Vc — Mj/^, the rj c width T/^ t , the mass resolution 
of the J/ip <Jj/^, the normalization and slope of the background and the number of events in 
the rj c and J/ip peaks. The result of the fit is presented in tabled 

Table 1: Result of the unbinned maximum likelihood fit of the mass spectrum. The resolutions of the J/ip and 

the r/ c peaks are o~j/^ and a Va , respectively . 





1 tot 

[MeV/c 2 ] 


Mj/tp - M Vc 
[MeV/c 2 ] 


[MeV/c 2 ] 


[MeV/c 2 ] 


[MeV/c 2 ] 


fit 


33.3 ±2.5 


113.6 ± 1.2 


3094.0 ± 0.8 


7.5 ±0.8 


oj/^ - 0.8 



1.3 Systematic errors 

A source of a systematic error on is the uncertainty of the detector mass resolution. This 
contribution is estimated to 0.4 MeV/c 2 by using the mass resolution for the rj c obtained from 
Monte Carlo events for the fit. The systematic error due to the parametrization of the back- 
ground is estimated by varying the mass range of the fit. A variation from 2.4-3.6 GeV/c 2 to 
2.7-3.3 GeV/c 2 changes the value of Tf ot by 0.7 MeV/c 2 . 

The reconstructed masses Mjj^ and M^ c are shifted in Monte Carlo simulation by -1.1 MeV/c 2 
compared to the values used as input. This bias does not affect the mass difference. However 
the fitted J/ip mass is still shifted by -1.8 MeV/c 2 from the well established PDG value^ of 
(3096. 87± 0.04) MeV/c 2 after correction of the bias observed in the Monte Carlo simulation. We 
take this shift as a systematic error on the r] c mass since we do not know its source. The unkown 
mass shift is not necessarily the same for the r] c and the J/ip because the angular distribution of 
their decay products are different. The shift could be caused by inhomogeneities in the magnetic 
field which may not be correctly modeled in Monte Carlo. 

1.4 Evidence for the r] c (2S) state 

The r/ c (2S) is not experimentally established yet^. Its theoretically predicted mass is in the 
range ^ 3583 to 3640 MeV/ c 2 . There exist two experimental results ^ ^ which are not in 
agreement. 

The upper right histogram of figurenshows the K®K if^ invariant mass spectrum extended 
to higher masses. An enhancement can be seen at about 3.63 GeV/c 2 which lies within the 
predicted mass region of the r] c (2S) and hence is interpreted as an evidence of this state. The 
mass spectrum is fitted using a Breit-Wigner function convoluted with a Gaussian to describe 
the enhancement. The standard deviation of the Gaussian is set to 8.5 MeV/c 2 which was 
estimated from a Monte Carlo sample. The fit results in T^ 2 "^ = (20 ± 10) MeV/c 2 and 
M Vc (2S) = (3633.3 ± 5.0) MeV/c 2 . The r] c (2S) peak contains 86 ± 23 events. 

The systematic error associated to the i] c (2S) width due to an uncertainty of the mass 
resolution is 1.6 MeV/c 2 . The background parametrization gives a contribution of 3.3 MeV/c 2 
to the systematic error on the total width. 
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Figure 1: Various histograms of the KgK ± TY^ invariant mass spectrum. The large peak around 2983 MeV/c 2 is 
the rjc, the peak next to it is the J/tp. An evidence of a i] c (2S) signal can be seen around 3630 MeV/c 2 . 

In summary, we have measured the mass and width of the rf c meson M Vc = (2983.3 ± 
1.2 (stat) ± 1.8 (syst)) MeV/c 2 and Yf ot = (33.3 ± 2.5 (stat) ± 0.8 (syst)) MeV/c 2 and observed 
an evidence of the rj c (2S). The mass and width of the r) c (2S) were measured to M Vc ( 2 s) = 

(3632.2 ± 5.0 ± 1.8) MeV/c 2 and T v t ^ 2S) = (20 ± 10 ± 4) MeV/c 2 , respectively. 

The measured value of the total width of the r\ c deviates from the world average by more 
than three standard deviations while the measured mass agrees well with the world average. 
The mass of the enhancement interpreted as the rj c (2S) falls in the theoretically predicted range 
of the r) c (2S) mass. 



2 Measurement of the r lepton lifetime 

A measurement of the r lepton lifetime, combined with the branching ratio of r into leptons, the 
r and \x mass and the [i lifetime, provides a test of lepton universality. This test is experimentally 
limited by the uncertainties of the r lifetime and its leptonic branching fraction. 
The mean r lifetime (r T ) is given by 

= ^ry^>> ( 2 ) 

where M T is the r mass, (P T ) the average r momentum in the center-of-mass frame and (A T ) the 
mean r decay length. The average momentum is calculated from a Monte Carlo sample. The 
mean decay length is measured using decays where one r decays into three charged particles 
while the other one decays into a single charged particle. 

The tracks of the 3-prong r candidate are projected onto the transverse plane, which is 
perpendicular to the boost direction of the e + e~ annihilation. The 3-prong vertex x^ p in the 
transverse plane is reconstructed and the transverse decay length X Tjt is calculated by: 

K,t = (X3p ~ X bs ) ■ P3 P (3) 

where x^s is the center of the luminous region and p^ p the unit vector of the 3-prong momentum 
in the transverse plane. This method of reconstructing the decay length was chosen because it 
is less dependent on uncertainties of the detector alignment. The decay length in space A r is 



calculated by dividing the transverse decay length by the sine of the polar angle of the 3-prong 
momentum in the center-of-mass system. Figure El shows the decay length distribution for r 
candidates. 

The mean decay length is calculated by averaging over all events. The averaging is done 
without any error weighting of the events. Instead an azimuthal weighting is done to achieve an 
uniform azimuthal event distribution. This method of averaging is more robust with regards to 
uncertainties of the position of the luminous region. 
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Figure 2: Measured decay length for 
data and Monte Carlo. 



The bias of this analysis method was obtained from a Monte Carlo sample. Its uncertain- 
ties are the main contribution (1.3 fs) to the systematic error of the r lifetime. The detector 
alignment is another major source of systematic uncertainties (0.7 fs). 

The result based on 30.2 fb' 1 is r T = (290.8 ± 1.5(stat.) ± l.6(syst.)) fs. It agrees well with 
the world average of (290.6 ± 1.1) fs. 
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